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Abstract 

A method to calculate the nuclear double beta decay {2uf3l3 and Qvf3(3) amplitudes within the 
continuum quasiparticle random phase approximation (cQRPA) is formulated. Calculations of 
the (ifi transition amplitudes within the cQRPA are performed for '''^Ge, ^'^"Mo and ^^''Te. A 
rather simple nuclear Hamiltonian consisting of a phenomenological mean field and a zero-range 
residual particle-hole and particle-particle interaction is used. The calculated 2vf3f3 amplitudes 
are almost unaffected when the single-particle continuum is taken into account, whereas we find 
a regular suppression of the Oi^/3/5 amplitudes that can be associated with additional ground-state 
correlations owing to collective states in the continuum. It is expected that inclusion of nucleon 
pairing in the single-particle continuum will somewhat compensate this suppression. 
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I. INTRODUCTION 



Neutrino oscillation experiments have proven that neutrinos are massive particles (see, 
e.g., Ref. [li]). However, the absolute scale of the neutrino masses cannot in principle be 
deduced from the observed oscillations. To determine the absolute neutrino masses down 
to the level of tens of meV, study of the neutrinoless double beta decay (Oi//3/3) becomes 
indispensable. Furthermore, this process, which violates the total lepton number by two 
units, is an experimentum crucis to reveal the Majorana nature of neutrinos 0, 0, 0, [y]- 

The next generation of experiments (GERDA, CUORE, SuperNEMO etc.) has a great 
discovery potential for observation of 0z//3/3 decay and for providing reliable measurements of 
the corresponding lifetimes. The determination of the effective Majorana mass (or relevant 
GUT and SUSY parameters depending on what mechanism of the Qv(3(3 decay dominates) 
from experimental data can be only as good as the knowledge of the nuclear matrix elements 
M^^ on which the {^v(5(3 decay rates depend. Thus, a better understanding of the nuclear 
structure effects important for describing the matrix elements is needed to interpret the data 
accurately. It is crucial in this connection to develop theoretical methods capable of reliably 
evaluating the nuclear matrix elements, and to realistically assess their uncertainties. 

At present, the most elaborate analysis of the uncertainties in the Qv(3(3 decay nuclear ma- 
trix elements M^^ calculated within the quasiparticle random phase approximation (QRPA) 
and the renormalized quasiparticle random phase approximation (RQRPA) has been per- 
formed in recent papers 0,0] (the bases comprising N =2, 3 and 5 major oscillator shells 
were used). The experimental 2i//9~/5~ decay rates have been used there to adjust the most 
relevant parameter, the strength gpp of the particle-particle interaction. The major observa- 
tion of Refs. [sl, 0] is that such a procedure makes the calculated M^'^ essentially independent 
of the size of the single-particle (s.p.) basis of the QRPA . Furthermore, the matrix elements 
have been demonstrated to also become rather stable with respect to the possible quenching 
of the axial vector coupling constant qa- 

The calculations in Refs. [g], 0] were performed within "the standard QRPA" scheme in 
which a discrete s.p. basis and the harmonic oscillator wave functions as the s.p. wave func- 
tions are used to build the BCS ground state and the spectrum of the excited states. Keeping 
in mind that many multipoles contribute appreciably to M^^ , one can a priori expect that 
enlargement of the model space should lead to more accurate matrix elements M^^ . (In 
other words, any basis truncation leads to an uncertainty.) This should be contrasted with 
the case of the 2vl3(3 amplitude to which only Gamow- Teller transitions contribute and a 
s.p. basis of N =1-2 major shells is good enough. In this respect, it would be interesting 
to test the stability of the calculated M^'^ found in Refs. [gI, 0] by letting ^ oo. Thus, if 
one could include the entire s.p. basis into the calculation scheme, the question about the 
dependence of the QRPA results on size of the s.p. basis as a source of the uncertainties in 
the calculated M^^ would become irrelevant. 

There is no problem within the QRPA for including low-lying major shells composed of 
bound s.p. states into the model space. But inclusion of major shells lying much higher than 
the Fermi level immediately encounters principal limitations of approximation of the contin- 
uum of unbound s.p. states by discrete levels. Basically, only one major shell, lying higher 
than the Fermi shell (already containing quasistationary states), can safely be considered. 

The only possible way to treat properly the s.p. continuum is provided within the con- 
tinuum quasiparticle random phase approximation (cQRPA). The continuum random phase 
approximation (cRPA) was formulated about 30 years ago in the pioneering work of Shlomo 
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and Bertsch |8| and since then has been used to successfully describe structure and decay 
properties of various giant resonances , muon capture fl^ , and neutrino-nucleus reactions 



with large momentum transfer [HI]. To apply the cRPA in open-shell nuclei one has to take 
nucleon pairing into consideration. This requires development of a quasiparticle version of 
the cRPA, namely, a continuum QRPA approach. Such a cQRPA approach to describing 



charge-exchange excitations has been developed in Refs. [12, 1131 . 

The cQRPA provides a regular way of using realistic wave functions of unbound s.p. 
states in terms of the s.p. Green's functions without the need to approximate them by 
the oscillator ones. Moreover, having an alternative formulation of the QRPA can help to 
understand current QRPA results and their deficiencies. 

Two principal effects of taking into account the s.p. continuum within the proton-neutron 
QRPA (pn-QRPA), which affect the calculated values of M^" in an opposite way, can be 
expected. First, additional ground-state (g.s.) correlations can appear because of collective 
multipole states in the continuum, which generally have a tendency to decrease M^". Second, 
pairing in the continuum can increase Oz//3/3 matrix element M^'^ (see the relevant discussion 
in Ref. Q). 

The principal aim of this work is to formulate for the first time a nuclear structure 
framework for calculating the double beta decay matrix elements M^^ and M^'^ within 
the cQRPA and to test within this method the stability of the calculated M'^^ found in 
Refs. p, 0|- As a first step, a simpler version of the cQRPA with nucleon pairing realized 
only on a discrete basis is applied in the present work; therefore, the calculated M°^'s of 
this paper should be considered lower limits for the matrix elements within the cQRPA. To 
consistently include nucleon pairing in the continuum within the cQRPA is a formidable 
task and is postponed to future publications. 

We merely focus here on a qualitative discussion of the relative effect obtained within 
the cQRPA in comparison with the standard discrete QRPA. Therefore, the M"*^ values 
of the present work may be somewhat different from those of Refs. [g], 0] since we have 
not implemented the most elaborate representation for the neutrino potential (modified 
by the finite nucleon size correction, higher order terms of the nucleon weak current, etc.; 
see, e.g., The nucleon-nucleon short-range correlations (SRC) are implemented here in 



the usual way, in terms of the Jastrow-like functions [15|. This, however, might lead to a 
overestimation of the effect of the SRC (see the recent discussion in Refs. HH). 

The paper is organized as follows: The pn-QRPA equations in the coordinate presentation 
and the way to take into account the s.p. continuum in them are given in the first two parts 
of Sec. II. In the latter two parts of that section formulas for calculating strength functions 
and M"^" and M^'^ are presented. In Sec. IIIII we present the results and we give conclusions 
in Sec. HVl 



II. CONTINUUM QRPA 

Since its formulation in the pioneering work of Shlomo and Bertsch [8|], the cRPA has 
long been used for to successfully describe structure and decay properties of various giant 
resonances and their high-lying overtones embedded in the single-particle continuum. The 
structure of the overtones is formed by the s.p. excitations changing the s.p. radial quantum 
number (which correspond to transitions over two or more major shells). Their contribution 
to the nuclear multipole response is marked if probe operators have nontrivial radial depen- 



dence, which is the case, for example, for muon capture [10| and neutrino-nucleus reactions 
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with large momentum transfer [ijj . The direct nucleon decay of various giant resonances and 
their overtones has been extensively analyzed within the cRPA by Urin and collaborators 
(see, e.g., Q). 

To apply the cRPA in open-shell nuclei one has to take nucleon pairing into considera- 
tion. This requires development of a quasiparticle version of the cRPA, namely, the cQRPA 
approach. The approach should account for the important influence of the residual particle- 
particle (p-p) interaction along with the particle-hole (p-h) one included within the usual 
cRPA. Such a cQRPA approach based on the coordinate space Hartree-Fock-Bogolyubov for- 
malism has been formulated and applied recently to describe strength functions of different 



multipole excitations without charge exchange [IT] . 

A pn-cQRPA approach to describing charge-exchange excitations was developed in 
Ref. and, independently, in Ref. In Ref. [l8| the approach was applied to analyze 
the low-energy part of the Gamow- Teller (GT) strength distribution relevant for descrip- 
tion of single beta decay in astrophysical applications. In [T^ the Fermi and GT strength 
distributions in semimagic nuclei were described within a wide excitation-energy interval 
that includes the overtones of the IAS and GTR, the so-called monopole and spin-monopole 
resonances. 

In describing the Oi>(3f3 decay, some transition strength into the s.p. continuum is missing 
within the standard QRPA calculation scheme, especially for the high-multipole excitations 
with L > 2 (compare, e.g., with the description of muon capture where the contribution of 
the highly excited giant resonances dominates The contribution of these multipoles 

to M^'^ becomes particularly important because the monopole (Fermi and Gamow- Teller) 
contributions are suppressed by symmetry constraints (see, e.g., the multipole decomposition 
of M^"" in Fig. 5 of Ref. 0; for a recent general discussion of how the SU(4)-symmetry 
violation by the residual p-p interaction affects M'^'^ see Ref. [191]). Thus, M^'^ gets strongly 
suppressed by the g.s. correlations, short-range correlations, etc.; therefore fine effects (such 
as influence of the s.p. continuum) can be expected to come into play. The first attempt to 
briefly describe the jS/S observables within the pn-cQRPA was undertaken in Ref. 13 . 



A. pn-QRPA equations in coordinate representation 

The system of homogeneous equations for the forward and backward amplitudes Xnu 
and Ynu , respectively, is usually solved to calculate the energies a;^ and the wave functions 



I J'^M, s) of excited states in isobaric odd-odd nuclei within the pn-QRPA (see, e.g., [2, l21 
here "s" labels the different QRPA states). However, it is impossible to handle an infinite 
number of amplitudes X, Y if one wants to include the continuum of unbound s.p. states. 
Instead, by going into the coordinate representation the pn-QRPA can be reformulated in 
equivalent terms of four-component radial transition density {qY^^^} (/ = 1, ... ,4) defined 
for each state | J'^M, s) . The components are determined by the standard QRPA amplitudes 
X and Y as follows: 
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r) = ^/?-X..(r), 



TDTTV 



(1) 



where ti 

,(JLS) 



coefficients of Bogolyubov transformation and 



and V are the 

^w'm with (vr) = {jj^) [(z/) = {jX)] and r-^Trl?^) ['^"^^.(j^)] being the s.p. 

proton (neutron) quantum numbers and radial wave functions, respectively. In Eq. ([1]) the 
spin-angular variables are already separated out since the nuclear response to a probe op- 
erator having definite spin-angular symmetry determined by the irreducible spin-angular 
tensor Tjlsm{^) is calculated, and = '^/^^{'^\\^JLs\\^) represents the corresponding 

reduced matrix element. Hereafter, we shall systematically omit the superscript ^^{JLS)" 
when it does not lead to a confusion. 

According to the definition of Eq. ([1]), the elements Qi, §2, Qs and can be called 
the particle-hole, hole-particle, particle-particle, and hole-hole components of the transition 
density, respectively, and can be generally considered as a four-dimensional vector: {f?/}. In 
particular, the transition matrix element to a state \s, JM) corresponding to a particle- hole 
operator 



yi-) . 



= V(jLS){ra)TjLSM{na) 



(2) 
(3) 



is determined by a one-dimensional integral of the product of the element Qi and the radial 
dependence of the operator Eqs. 



(-) 



JLSM 



|0) 



g[^^^\r)V(^jLS){r) dr. 



(4) 



The pn-QRPA system of integral equations for the elements Q^"^^^^ follows from the stan- 
dard pn-QRPA equations for the X and Y amplitudes (see, e.g., Refs. 0, iH) by making 
use of the definition from Eq. ([1]): 



{JLS)f 



K 



(5) 



Here, 



(rr') ^F^^^\r, r') is the radial part of the residual interaction in the K channel (where 
K = 1,2 for the p-h channel, ii' = 3,4 for the p-p channel and the so-called symmetric 



approximation is used here). The 4x4 matrix 



AiK{rir2,uj) is the radial part of the 
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free two-quasiparticle propagator (response function): 



AiK{nr2^uj) = E AYK{uj)x.u{n)x.u{r2)-. AY, = AYk (6) 
ul vl \ ( vl ul 



'-^ — E^^y —io — Ej^i, J \uj — Et,j, —10 — Et,^ 

^12 — ^^34 — ^ r 



LU — Ej^y —UJ — Ej^y 

A-,^{uj) = AY{-oo), A-H = A-(-.;), Ar,{uj) = AY{-oo), = 

with i^TT!/ = E-,, + Ey, where E.„ and Ei, are the proton and neutron quasiparticle energy, 
respectively. The expressions for the elements of the free two-quasiparticle propagator Ajk 
can also be obtained by making use of the regular and anomalous s.p. Green's functions for 



Fermi systems with nucleon pairing (see, e.g., [13|), in an analogous way to that described 
in the monograph [2^ for response of Fermi systems to a s.p. probe operator acting in the 
neutral channel. 

These equations allow a compact schematic representation when the spin-angular vari- 
ables are not separated out. In such a case, the substitutions — > x^, Xa{r) Xa(x), 
g ^ Q, Ajk — > A/ii-, and Fk — > have to be made and the factor t^'^^^^i should be 
omitted in the formulas. Then, schematically denoting the double integration over x', x" in 
eq. ([5]) as {. . . }, one can rewrite the equation as 

Qj = {AikFkQk}, (7) 

where summation over the repeated index K on the right-hand side is assumed. 

The total two-quasiparticle propagator (two-quasiparticle Green function) A that includes 
the QRPA iterations of the p-h and p-p interactions is very useful in practical applications. 
It satisfies an integral equation of the Bethe-Salpeter type Ajk = Ajk + {Ajk'Fk'Ak'k}'- 

if/^)(rr', u) = AY^^'\rr', u;) + E / u;) 4''^''Vi^2) A^K^K\r2r', u) dndr,.{^) 

K' 

The following spectral decompositions hold for the radial components Aii{uj) , Ai2{uj) 
and A22{'^y 

in(nr„ c.) = y MiiiMM _ y M(!iMM 

s s 

A22irir2, uj) = An{rir2, -u) (9) 
4 ^ Qliri)02ir2) ^2(^)^1 (^2) 

^-^ UJ — UUg + to UJ + UJs — 10 

s s 

These components are the only ones that will be needed in the following for description 
of single and double beta decay transition probabilities. The spectral decompositions for 
the other elements of A can readily be written down in analogy to Eqs. (Q. Thus, one sees 
from Eqs. ([9]) that all the information about the QRPA solutions [energies ujg and transition 
densities ^?7,s(^)] resides in the poles of the total two-quasiparticle propagator A. 
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In this paper, the residual isovector particle-hole interaction and the particle-particle 
interaction in both the neutral (pairing) and charge-exchange channels are chosen in the 
form of the Landau-Migdal forces of zero range (proportional to the spatial 6 function) j22| , 



which is similar to the choice of Refs. |23l. |24|| . The effective isovector particle- hole interaction 
Fk (for K = 1,2) is given by 

F^(xi, X2) = Co(/;, + /ifii ■ 62)ri ■ T,6{r, - r^), (10) 

where /^^^ and /^^^ are the phenomenological Landau-Migdal parameters. Hereafter, all the 
strength parameters of the residual interactions are given in units of Cq = 300 MeV- fm^. 
The residual p-p interaction (for K = 3,4) is given by a similar expression: 

F^(xi, X2) = -Coig^p + glpQ, ■ Q2)S{r, - r,), (11) 

and the pairing interaction: 

FP'^^^(xi, X2) = -Cog^'^'^-Sir, - r2). (12) 

The pairing strengths g^"'^'^ and g^""^"^ for neutron and proton subsystems are fixed within 
the BCS model to reproduce the experimental neutron and proton pairing energies. All 
the other strength parameters in the particle-particle channel are always given relative to 

{gr'+grn/^- 



B. Taking the single-particle continuum into consideration 

The coordinate-space version of the pn-QRPA described in the preceding section is es- 
pecially suitable for taking the s.p. continuum into consideration. But before proceeding 
with the continuum, it is worth noting that if one lets the double sums in Eq. run just 
over finite sets of proton and neutron s.p. states, the presented version of the pn-QRPA 
is fully equivalent to the usual "discrete", one, which is formulated in terms of X and Y 
amplitudes. We make use of this fact in order to check the calculation scheme by comparing 
"discrete" QRPA results calculated in these two different, but formally equivalent, ways. As 
anticipated, the results are the same within the accuracy of the numerical techniques used. 

To take the s.p. continuum into consideration, the double-sum representation for the free 
response function should be transformed according to the following prescription: 

1. The Bogolyubov coefficients Va, Ua and the quasiparticle energies Ea are approximated 
by their no-pairing values Va = 0, Mq = 1, and Ea = Sa — Aj for those s.p. states in the 
s.p. continuum that lie far up of the chemical potential Aj [i.e. — Aj ^ A^, "i" = 
"p" (protons) or "n" (neutrons)]. The accuracy of this approximation is "^^ich 
is good enough already for — Aj > E^ax — 10 MeV. The usual BCS representations 
for Va, Ua and Ea are taken for all the other s.p. states with Ea — \i < Emax- 

2. The radial single-particle Green's function 5'(Q)(rir2, e) = Xa( i)Xa{ 2) ^gg^ 

explicitly perform the sum over the s.p. states in the continuum. Here, the sum 
runs over different radial quantum numbers for a given spin-angular symmetry (a). 
The Green's function satisfies the inhomogeneous radial s.p. Schrodinger equation 
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[^o(a)(^) " £]9ia){rr' , e) = —S{r — r') and can be constructed as a product of regular 
and irregular solutions of the homogeneous equation [/io(a)('") — 
(see, e.g., Refs. BiH). 



reg,irreg , 



As a result, we get from Eq. ([6]) the following representation for the components Ajk of 
the free two-quasiparticle propagator: 



X^v{,ri)xnu{,r2) + ^ Xu{.n)Xu{,r2) 9[n){nr2, Xp + uj - E^) 

(13) 



Ai2(rir2,cj) 



uj - E^ 



+ 



1 



-uj - E^ 



i/<,7r< 



uz 



u — Et, 



+ 



-cij — E^ 



Xivu{ri)x^i,{r2) 



i^<,(71-) 



^33 (^^^2,^^) = 



UJ - E^ 



+ 



-cj - E^ 



Xttv 



!^<.(t) 

+ X (^Wh) M^X^(n)X7r(r2)^(^)(rir2,A„ + w -E^) 



7r<,('^) 

v444(u;) = As3{-uj); Au{uj) = Ais{-uj,7r ^ u); ^24(0;) = ^i3(-t^); v423(c^) = ^13(0;, tt ^ z/) 



where 7r< (z/<) means vr < Umax < ^'max), where Umax {^max) is the s.p. state with the 
largest energy included into the BCS basis for which E^^^^ = e^^^^-Xp {E^^^^ = e^^^^-Xn) 
with the acceptable accuracy as previously described, and 



g[„)irir2,e) = g(^)inr2,e) - ^ 



X7r(ri)x^(r2) 



(14) 



is the subtracted radial s.p. Green's function (the Green's function from which the contri- 
bution of all discrete s.p. states and those quasidiscrete s.p. states included in the BCS 
basis is subtracted). 



C. Strength functions 

Different strength functions can be readily calculated in terms of the imaginary part 
of the total two-quasiparticle propagator ImA. The strength function corresponding to a 
charge-exchange single-particle operator V^jj^ acting in the [3^^'^ channel, 

VjVsM = J2"^JLSMMT!,'^\ (15) 
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where y^jLSAii^a) is given by Eq. ([3]), is defined by the usual expression: 

s 

with ujf = Ef — Eq being the excitation energy of the corresponding isobaric nucleus 
(A^ =F 1, Z ± 1) relative to the ground state |0) of the parent nucleus {N, Z) with energy Eq. 
Making use of the spectral decomposition ^ one can easily verify the following integral 
representations of the strength functions: 

S^-\uj-) = -i jvAnMv} = -ilm/V(ji5)(ri)iSf''^(rir2;cu)V(ji5)(r2)rfrif/r2,(16) 
= -i |VA22(^)V} = -ilm/ V(jz.5)(ri)4f ''^(rir2;^)V(ji5)(r2)rfrirfr2.(17) 

where uo^ = ± (Ap — A„) ± (m„ — m^) represents the calculated excitation energy relative 
to the g.s. of the parent nucleus and rrip, m„ are proton and neutron masses. The pn-QRPA 
excitation spectrum, originally calculated in terms of a;, gets a constant energy shift to be 
represented in terms of cj^, because the modified nuclear Hamiltonian H — XpZ — XnN (as 
in the BCS model) is used within the QRPA and the model nuclear Hamiltonian does not 
contain the rest energies of nucleons. 

One can also define a nondiagonal strength function as 

^-^M = 5^(Olvj^)^|J-M,.)(J-M,.|Ki-^^|0)<5(^-a;:) (18) 

s 

with Us = Eg — {Eq + Eqi)/2 = {uj~ + uj'^)/2 and Vj^gj^ being the time reverse of VjlIm- 
Such a strength function is closely related to the amplitude of the 2u(3(3 decay, when |0) 
and |0') are the g.s. wave functions of the initial (decaying) and final (product) nuclei, 
respectively. 

To calculate Sy \uj) within the pn-QRPA one faces the usual problem that the spectrum 
\s) comes out slightly different when calculated with respect to |0) or |0'). Identifying the 
QRPA vacuum |0') with that of |0), one gets Us = ojs and 

S^-\uj) = -- {vAi2(a;)v} = --Im / V(jL5)(ri)4f''^(rir2;cu)V(jL5)(r2)f/rirfr2 (19) 
or, alternatively, identifying |0) with |0') 

S^-\uj) = |vA;2(a;)v} = -l^^j V^JLS){ri)A''-/^'''\nr2-,u)V^jLs){r2)dndr2 (20) 
where A' is calculated with respect to the g.s. |0') of the final nucleus. 



D. Description of (3(3 decay within the cQRPA 

The spectral decomposition of the two-quasiparticle propagator A Qj can be used for 
calculation of PP decay matrix elements in a similar way as described in Sec. Ill Cl for \uj) 

(Hi. 
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The 2uPP decay amplitude Mq^ is defined by the following expression: 



= V ^^"^ "zr"^ (21) 



(0'||G(-)||s)(s||G(-)||0) 



(jj 



where G^-^ = ^^B„ r" and again cj, = - (£"0 + Eo')/2 = {uj; + uj'^)/2. 

By using the spectral decomposition of Eqs. ([9]) for {JLS) = (101) and the approximation 
that the QRPA vacuum |0') of the final g.s. is the same as |0) of the initial g.s. (the same 
approximation as used in Refs. [2^, 24|), the amplitude fl2Tl) is simply given by one-half of 
the corresponding static nuclear polarizability with respect to the external s.p. field fir^: 

M^^^ = -i |6Ai2(cj = 0)fi} = -Gtt J A^^^^^\nr2]uj = 0)dndr2 

where lJs = ujg is used consistently in this approximation. 

The same procedure can be applied to calculate within the cQRPA the matrix element 
(0'|iy*^ -"lO) of a two-body scalar operator 

W.(--) = ^W(x.,x,)ri-)r(-) (22) 

ab 

W(x„,X6)= W(^jLs){ra,rb)TjLSM{na)T}Ls^j{nb) 



JLSM 



between the ground states |0) and |0'). It is given by a sum of all partial contributions 

M(— ) = -i / c/u; Im | WAi2(u;) } = Y1 ^^'^^^^ (^3) 

JL 

MiJLS) ^ j du j iy(ji.5)(ri,r2)Im4f ^)(rir2;o;)ciridr2 (24) 

(where the identification of the ground states as described previously has to be done). 

The neutrino potential appearing in the description of the Oz//3/3 decay (see,e.g., Refs. [2I, 
i): 

„h V 9 a/ 



in the simplest Coulomb approximation Py{\ra — ri,\ = Tab) = R/^ab (with R = 1.23 A^/^ 
fm being the nuclear radius) has the well-known partial radial components W(^jLS){ra, fb) = 

2lTT?> ) (^^^'^^ ^< ~ Kiin(T'a, rb), r> = max(ra, r^)). When one takes into account both 
the energy dependence of the neutrino potential Pi, = P^{rab,uj) and the usual modification 
P„ Pvfj with the Jastrow-like function fj{rab) = 1 — e^'^^'"^!'(l — ^2r1b) to account for the 
SRC of the two initial neutrons and two final protons , the decomposition of the neutrino 
potential over the Legendre polynomials P^ can be done numerically: 

TT 

W"(jLS)(^b^2,cu) = j dcoseuPL{cosei2)Pu{ri2,uj). (26) 
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In this first application of tlie cQRPA tlie corrections from the high-order terms in the 
nucleonic weak current and the finite nucleon size are omitted, which can lead to a slight 
overestimation of the calculated M^" by 20-30 % 0- 

Note that, within the cQRPA, in contrast to the discrete QRPA, one does not get an 
explicit set of QRPA energies and the energy integrations in the expressions for (I23p . (|2^ 
have to be performed on a mesh. For each point in the energy mesh the cQRPA equations 
dH]) are solved by discretizing the spatial integrals, thereby transforming them to a matrix 
representation. All this makes the calculation of M'^'^ and rather time consuming. 
Implementation of adaptive integration methods helps to optimize the integration over the 
energy. 



III. CALCULATION RESULTS 

We perform the first calculations of the jSjS transition amplitudes M^*^ and M^'^ within 
the cQRPA for ™Ge, ^°°Mo and ^^°Te. We also compare the results obtained with those 
calculated within the usual "discrete" version of the QRPA to see the influence of the single- 
particle continuum. 

For the first calculations of M^*^ and M^" within the continuum QRPA we adopt a rather 



simple nuclear Hamiltonian similar to that used in Refs. [23|, |2J] . The chosen nuclear mean 
field U{x) consists of the phenomenological isoscalar part Uo{x) along with the isovector 
Ui{x) and the Coulomb Uc{x) parts, both calculated consistently in the Hartree approx- 
imation (see Ref. The residual isovector particle- hole interaction and the particle- 
particle interaction in both the neutral (pairing) and charge-exchange channels are chosen 
in the form of zero-range forces [Eqs (fT0|) - [T2|) ]. All the strength parameters of the residual 
interactions are given in units of 300 MeV-fm^. 

The results calculated within the discrete-QRPA, labeled A, B, C, refer to different s.p. 
bases used. Case A corresponds to the large s.p. basis in the calculations: 16 successive 
s.p. levels comprising = 1 — 5 major Saxon- Woods shells for ^^Ge and ^°°Mo and 22 
successive s.p. levels (all bound s.p. states for neutrons and all bound s.p. states along with 
6 quasistationary states for protons) comprising = 1 — 6 major shells for ^^°Te. Note 
that the same s.p. basis is used within the cQRPA as the basis of the BCS problem. Case 
B corresponds to the small s.p. basis and is obtained from A by subtracting the six lowest 
s.p. levels comprising = 1 — 3 major shells (inert core of ^°Ca). Case C corresponds 
to the solution of the QRPA equations in the large basis A, in which, however, the BCS 
problem is solved in the small basis B (i.e. the Bogolyubov coefficients f „ = 1 are taken for 
the six lowest s.p. levels). This tests the approximations involved in the current cQRPA 
calculations, namely, that Va is set exactly to zero for the s.p. states lying above the BCS 
basis. 

Fixing the model parameters is done in the following way: 

• The p-h isovector strength /^^ is chosen equal to unity, /^^ = 1.0. This allows repro- 
duction of the experimental nucleon binding energies for closed-shell nuclei by imple- 
menting the isospin self-consistency of the symmetry potential Ui{x) of the mean field 
with the isovector p-h interaction (see Ref. 13|). 



The p-h spin-isovector strength /^^ is fitted to reproduce the experimental energy of 
the GTR. 
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TABLE I: Choice of the model parameters. See text for explanation of the choices A, B, C for 
different model spaces within the discrete QRPA. The two lines for the parameters of for each 
PP decay are fitted by using the QRPA set of the intermediate states built with respect to the g.s 
of the initial nucleus (e.g. ''^Ge) and the g.s. of the final nucleus (e.g. ™Se). The p-h strength 
fpp is adjusted to reproduce the energy of the Gamow- Teller resonance in the initial nucleus (e.g. 
76Ge). 



Nuclear 




strength 


discrete QRPA 


continuuni 


transition 


MeV-1 


parameters 


A 


B 


C 


QRPA 






9pp 


0.98 


0.97 


0.97 


0.91 


76Ge^76se 


0.15 


9pp 


1.10 


1.10 


1.10 


1.01 






fph 


0.40 


0.43 


0.43 


0.40 






alp 


1.28 


1.31 


1.30 


1.10 


ioomo^iooRu 


0.24 


9pp 


1.43 


1.50 


1.49 


1.23 






fph 


0.70 


0.70 


0.70 


0.75 






9pp 


1.23 


1.25 


1.25 


1.11 


i30Tg^i30xe 


0.03 


9lp 


1.25 


1.28 


1.28 


1.13 






fph 


0.60 


0.60 


0.60 


0.63 



• The pairing strengths (7^°*^, Qp"''^'^ are fixed within the BCS model to reproduce the 
experimental pairing energies. As aheady mentioned, all the other strength parameters 
in the particle-particle channel are given relative to (5'^"*'^ + g^°'^^)/2. 

• By choosing the p-p isovector strength Qpp = 1.0 we restore approximately the isospin 
self-consistency of the total residual p-p interaction. 

• The p-p spin-isovector strength gL^ is chosen to reproduce the experimental (positive) 
value of M^'^ as done in Refs. 0,0]. 

In Fig. 1 the calculated g^p dependence of M"^^ is plotted. Calculated M"^^ depends on 
the choice of the g.s. (either initial or final) with respect to which the QRPA equations are 
solved. It can be seen in the figure that the M'^^ calculated within the discrete QRPA and 
the cQRPA are almost the same at small values of g^p, though the difference becomes visible 
while gpp grows. Thus, the correction to the M^^ coming from the s.p. continuum is small, 
as can be expected for the Gamow- Teller transitions. 

The values of the strength parameters /^^^ and g^p fixed according to our prescription 
are listed in Table [I] for both the discrete and continuum versions of the QRPA. Because 
calculated M'^^ depends on the choice of the g.s. (either initial or final) with respect to which 
the QRPA equations are solved, two sets of g^p are obtained. The upper and lower lines 
for each decay sequence in Table [T] contain g^p fitted for initial and final nuclei, respectively. 
One sees that the difference in the obtained g^p is almost negligible for ^^°Te— ^^^''Xe, but it 
becomes 6gl^ ^ 0.2 for i°°Mo^^°°Ru. 

The calculated values of M*^^ are given in Table [TTl for both versions of the QRPA for 
gA = 1.25. The numbers in parentheses are the matrix elements calculated with inclusion 
of the SRC in terms of the Jastrow-like function. The two lines of results for each (5(3 decay 
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FIG. 1: (Color online) The calculated dependence of the 2v(3(3 matrix element on the p-p strength 
Qpp in both discrete (set A, dashed lines) and continuum QRPA (solid lines). The calculations 
performed by using the QRPA set of the intermediate states built with respect to the g.s of the 
initial nucleus(e.g. ''^Ge) and the g.s. of the final nucleus (e.g. ^^Se) are depicted by thick and thin 
lines, respectively. The solid horizontal line gives the (positive) experimental value of the 2uf5(5 
decay matrix elements used in the calculation to fit gl^ . 



chain contain M^^ calculated with respect to the g.s. of the initial or final nucleus in the 
decay (see Sec. Ill D| ). In the present first application of the cQRPA, neither the finite nucleon 
size nor the higher order terms of the nucleon current are considered. (They usually bring 
an additional reduction of M^" by about 30%; see, e.g., Ref. 01). 

The contribution of the multipoles with L = — 5 are included in the calculations of M^" . 
The contributions with L > 5 (which can increase M^" in total by about 10%) are omitted 
here as the corresponding parts in the transition operator probe the short-range behavior 
of the nucleon-nucleon wave function that cannot be well described within the QRPA. It is 
known that the RPA in medium is formulated to describe propagation of small-amplitude 
density fluctuations and only the ring diagrams are summed (see, e.g., Ref. [25|). This is 
a quite suitable approximation to deal with collective long-wave excitations, but for the 



13 



short-range ones the diagrams that are left out of the RPA method become important. 

As one sees from comparison of the discrete QRPA results listed in columns A and B, the 
calculated M'^'^ values for different basis sizes come out very close to each other provided 
the p-p interaction parameter g^p is fixed to reproduce the experimental 2z//3/3 decay matrix 
elements Mf^^. This result provides an independent confirmation of the main conclusion of 
Ref. 0] which is obtained here for a different nuclear Hamiltonian and by solving the pn- 
QRPA equations in the coordinate representation. Also, the M^'^ values obtained by using 
only initial QRPA g.s. (upper line) or final one (lower line) in the calculation are rather 
closer to each other. 

The matrix elements M^'^ calculated by taking into account the SRC in terms of the 
Jastrow functions get suppressed by about 25-30 % in the present calculation, that is in 
quantitative agreement with other recent calculations 0, Eil but do not support the old 



results of Ref. [24| where the strong suppression was found. However, one should keep in 



mind that this way of describing the SRC can be rather rough and can lead to overestimation 
of the suppression of the M^'^ . Other methods of describing the SRC such as UCOM [2^ 
give much less suppression in the calculated M (only about 10 %) and this important issue 
is currently under intensive study [13, 3 • 

The matrix elements M'^^ calculated within the cQRPA (the last column of Table HTll . 
they are systematically smaller than the discrete QRPA ones (columns A and B). The 
suppression varies from about 30 % for ^^Ge to a factor 2 for ^"^Mo and ^^°Te. The origin 
of this suppression can be associated with additional ground state correlations appearing 
because of highly excited collective states embedded in the s.p. continuum. Transitions to 
these states are naturally described within the cQRPA in terms of the s.p. Green's functions 
(see Sec. Ill Bp . However, the applied version of the cQRPA in this work does not include 
nucleon pairing in the s.p. continuum and therefore the values obtained here should 
be treated only as lower limits. Inclusion of nucleon pairing in the s.p. continuum (which 
is a formidable task) will definitely lead to an increase of the matrix elements within the 
cQRPA. 

To demonstrate the importance of nucleon pairing far from the Fermi level for quantitative 
description of the M^'^, the numbers listed in the column C of Table [III can be compared 
with those in columns A and B. Case C is introduced, as previously described, to test within 
the discrete-QRPA the neglect of pairing far from the Fermi level, in a manner similar to 
how it is done in the present version of the cQRPA. Inspecting column C, one sees a marked 
reduction, by about 30 %, in the calculated Ou^P matrix elements. Thus, expanding the 
discrete QRPA basis from the "small" one of B to a "large" one of C which neglects pairing 
effects in the inert core but allows transitions from the inert core, leads to a suppression in 
the M^^ because of more g.s. correlations. The suppression, however, gets almost completely 
compensated as nucleon pairing is switched on in the inert core and one goes from case C 
to case A. The same sort of compensation is natural to expect in the case of the cQRPA 
when nucleon pairing is switched on in the single-particle continuum. However, one cannot 
exclude that the compensation is incomplete. 

Let us conclude with some words about possible prospects for taking nucleon pairing in 
the s.p. continuum into consideration within the approach described in here. Though possi- 
ble ways of treating the continuum pairing within the QRPA can be found in the literature 
(see, e.g., Ref. [17|), direct implementation of them would drastically increase the corre- 
sponding calculation efforts. One would need first to calculate the solutions u{r) and v{r) 
of the coordinate Hartree-Fock-Bogolyubov equation for positive energies, from which then 
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TABLE II: {)v[3[3 nuclear matrix elements evaluated without and with (in parentheses) the SRC in 
both discrete and continuum QRPA. See text for explanation of the choices A, B, C for different 
model spaces within the discrete QRPA. The two lines of results for each /?/? decay contain M^" 
calculated by using the QRPA set of the intermediate states built with respect to the g.s of the 
initial nucleus (e.g. ''^Ge) and the g.s. of the final nucleus (e.g. ^^Se). 



Nuclear 
transition 


discrete QRPA 


continuum 
QRPA 


A 


B 


C 


76Ge 


5.95 (4.30) 
5.44 (3.86) 


5.63 (4.19) 
5.22 (3.82) 


4.30 (3.19) 
3.81 (2.76) 


4.30 (3.09) 
3.63 (2.46) 


iooMo 


5.52 (3.88) 
4.19 (2.73) 


5.35 (3.84) 
4.00 (2.65) 


4.24 (3.00) 
2.91 (1.84) 


2.49 (1.67) 
1.39 (0.67) 


^ i30xe 


3.17 (2.19) 
4.69 (3.21) 


3.14 (2.20) 
4.67 (3.21) 


2.56 (1.78) 
3.77 (2.57) 


1.70 (1.12) 
2.03 (1.28) 



additional continuum contributions to the expressions (fT3|) for the response function should 
be constructed by direct integration over energy. Probably, a more economical way is to 
discretize the continuum by putting the nucleus in a large box. These further developments 
are postponed to a future publication that should then finally answer the question about 
stability of M^'^ with respect to the basis size. 

IV. CONCLUSIONS 

A continuum-QRPA approach to calculation of the nuclear double beta decay 2z//5/3- and 
Oi//3/3-amplitudes has been formulated. Calculations of the amplitudes M'^^ and M'^'^ within 
the cQRPA are performed for ''^Ge, ^°°Mo and ^^"Te. A rather simple nuclear Hamilto- 
nian consisting of phenomenological mean field and zero-range residual particle-hole and 
particle-particle interaction is used. The M"^" are almost unaffected when the single-particle 
continuum is taken into account. In contrast, we find a regular suppression of the Oi//5/? 
amplitude that can be associated with additional ground state correlations owing to collec- 
tive states in the continuum. The calculated M^'^ values of this paper should be considered 
as lower limits for the matrix elements within the cQRPA, as nucleon pairing is realized 
only on a discrete basis within the present version of the cQRPA. It is expected that future 
inclusion of nucleon pairing in the single-particle continuum will somewhat compensate the 
observed suppression of M^^ values. 
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